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■ Abstract 

In this paper we focus on the finite n probabihty distribution function of the largest 
eigenvalue in the classical Gaussian Ensemble of n x n matrices (GE„). We derive the 
finite n largest eigenvalue probability distribution function for the Gaussian Orthogonal and 
Symplectic Ensembles and also prove an Edgeworth type Theorem for the largest eigenvalue 

■ probability distribution function of Gaussian Symplectic Ensemble. The correction terms to 
the limiting probability distribution are expressed in terms of the same Painleve II functions 
appearing in the Tracy- Widom distribution. 



1 Introduction 

In applications of the limiting probability distributions laws from Random Matrix 
Theory (e.g. [1], [15]) it is important to have an estimate on the convergence rates, if 
possible have a control on this rate of convergence. For recent reviews of applications 
of these distributions we refer the reader to [HI d [HI [HI [25]. In our desire to control 
the rate of convergence of the probability distribution of the largest eigenvalue from 
the Gaussian Orthogonal Ensemble GOEn, Gaussian Unitary Ensemble GUEn and 
Gaussian Symplectic Ensembles OS En, we introduce a fine tuning constant c in the 
scaling of the desired eigenvalue. The finite large n expansion is therefore a function 
of c. We use this constant to fine-tune the convergence rate. This work was done for 
the GUEn and GOEn- In completing the same work for the GSEn, we decided to 
find a closed formula (opposed to a large n asymptotic formula) whenever possible 
for each function appearing in the the final expression for these probabilities distri- 
bution functions. This approach has the advantages that the final results is a finite 
n representation of the distributions functions, we only need to perform the large n 
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expansion once for each function to have an Edgeworth type expansion for the desired 
probabihty. 

The Gaussian /3-ensembles are probabihty spaces on n-tuples of random variables 
{Ai, . . . , An} (think of them as eigenvalues of a randomly chosen matrix from the 
ensemble.) with probability density that the variables lie in an infinitesimal intervals 
about the points Xi, . . . ,Xn is 

P„^(xi, ■ ■ ■ ,Xn) = C„^exp I J JJ |xj - Xfcl^, (1.1) 

\ 1 / j<k 

with 

— oo < < oo, for z = 1, ■ ■ ■ , n. (1.2) 

Here C„,^ is the normalizing constant such that the total integral over the Xj's is one. 
When /3 = 1 we have the GOE„, when /3 = 2 we have the GUE„ and when /3 = 4 we 
have the GSE„. We denote the largest eigenvalue by X^ax^ ^y 

F„„p{t) = ¥{Xi^ < t) (1.3) 

his probability distribution function. 

When /3 = 2, the harmonic oscillator wave functions (see [13], or [18] for a complete 
definition) 

V^.(a:) = ^^^^y^=y^i/.(x)e-V2 = 0, 1, 2, . . . 
play an important role. We also have the Hermite kernel 

Kn,2[x,y) = y_^^k[x)^k[y) = \ - , 1-4 

^-^ V 2 X — y 

which is the kernel of the integral operator Kn^2 acting on L^(t, oo) , with resolvent 
kernel 

Rn,2ix,y;t) = {I - Kn,2)'^ ■ Kn,2{x,y). (1.5) 

Note here that the dot denotes operator multiplication. We have the following rep- 
resentation of (11.11) . (see for example, [17] for a derivation of this result) 

P„,2(Xl,--- ,Xn) = det{Kn,2{Xi,Xj))i<ij <n. 

Following Tracy and Widom in[20l[2ll[22l[23l[2ll[25l[26],we define 

1 1 

(77,\ 4 f n \ 

2] ^n{x), i){x) = l-j (/5„_l(x), (1.6) 

by e the integral operator with kernel 

et{x) = ^sgn(x-t), (1.7) 
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D the differentiation with respect to the independent variable, 

g„„,(x;t) = ((J-i^„,2)-\xV) (1.8) 

and 

P„,,(x;t) = ((/-K„)-\a:V). (1.9) 
We introduce the following quantities 

QnAt) = Qn,^(t;t), Pn,^(t) = Pn,i(t; t) (1-10) 

UnA't) = {Qn,i,f), Vn,i{t) = {Pn,i,f), (l.U) 
Vn,i{t) = {Qn,i,^), and Wn,i{t) = {Pn,i,^). (1.12) 

Here (■,■) denotes the inner product on L^(t, oo). In our notation, the subscript 
without the n represents the scaled limit of that quantity when n goes to infinity, and 
we dropped the second subscript i when it is zero. 

If Ai is the Airy function, the kernel Kn,2{x, y) then scaled] to the Airy kernel 

Our conventions cirG clS follows: 

Q,ix;s) = {{I-KAir\x'Ai), Qoix;s) = Q(x;s), (1.14) 

Pi{x;s) = ((/-AaO^\ x^Ai'), Poix;s) = P(x;s), (1.15) 

Qiis) = Qi{s;s), qo{s) = q{s), pi{s) = Pi(s;s), pq{s) = p{s), (1.16) 

Ui{s) = {Qi,Ai), uo{s) = u{s), Vi{s) = {Pi, Ai), vo{s) = v{s), (1.17) 

Vi{s) = {Qi,Ai), vo{s)=v{s), Wi{t) = {Pi, Ai'), and wo{s)=w{s). (1.18) 

Here ( ■ , ■ ) denotes the inner product on L^(s, oo) and z = 0, 1, 2, ■ ■ ■ 

We also note that q{s) is the Haskins-Macleod solution to the Paileve II equation 

q"{s) = sq{s) + 2q^{s) with the boundary condition q{s) ~ Ai(s) as s — )■ oo. 

We use the subscript n for unsealed quantities only. 



t 

Qn{x]t)dx, 



7^„,l := / Rn{x,t;t)dx, Vn,i ■= / Pn{x;t)dx, Qn,l ■= / 

J — oo J — oo J ~ 

(1.19) 

and 

/oo poo 
et{x)Rn{x,t;t)dx, Vn,A{t) := et{x)Pn{x;t)dx, 
oo J — oo 

/oo 
et{x)Q„{x;t)dx. (1.20) 
-oo 



-'^as n oo in the change x ~ ■\/2(n + c) + 2 e X and y = ^^ (n + c) + 2 an bY 
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The epsilon quantities are 

Qn,e{x;t) = {{I -Kny\x,y),eip{y)), = g„,.(t;t) (1.21) 



UnA't) = (<5n,£(a;;i), V5(a;)), Vn,e{t) = {Qn,e{x] t) , ^J^x)) , (1.22) 

where ( ■ , ■ ) denotes the inner product on L^(t, oo). 

The GOE„ and GSE„ analogue of fll.27p in Theorem 11.11 bellow will follow from 
representations (equations (40) and (41) of [22].) 

FnAtf = Fn,2{t) ■ (^{l - Vn,e{t)) (l " ^^n,lW) " ^(w(t) " C^)Pn,l(t)) (1-23) 

and 



Fn,4{t/V2Y = Fn,2{t) ■ {{l - Vn,e{t)) (l + ^7^„,4(t)) + ^g„,.(t) Pn,4(t)) • (1-24) 

To complete the work, we will use the close form of 7?.„,i, Vn,i, 'R-n,^, 'Pn,^, Vn,e, andg^.e 
derived in [5] to give a finite n representation of and -Fn,4, then the corresponding 
large n expansion to find an Edgeworth type expansion for Fn,4 as outlined in |4]. We 
will also need the following results, 

Theorem 1.1. f3] If we set 

r(s) = (2(n + c))^ + 2~^n^^ s and (1.25) 

Ec,2 '■= Ec,2is) = 2wi - 3u2 + (-20c^ + 3)vo + UiVq - UoVi + UqVq - ulwo- (1-26) 
Then as n ^ oo 

Fn,2{r{s)) = F2{s) |l + cuo{s) 72^3 _ ^^c,2(s) u-'i^ + 0{n-') (1.27) 

uniformly in s, and 

F2(s) = Jimi^„,2(t) = ^^p(-^ {x-s)qixYdx^ (1.28) 

is the Tracy- Widom distribution. 

To state the next result we need the following definitions 



a := a{s) = / q{x)u{x)dx, (1-29) 

J s 

/oo 
q{x)dx, (1.30) 
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p{x)dx = a{s) — q{s), 



(1.31) 



r] := r/(s) 



20^2 



(6qv + 3pu + 2p2 + 2piv + 2pvi — 2q2U — 2qiUi — 2qu2) [x] dx 



20c^q'{s) + 3p{s) 
20V2 



(1.32) 



EcAs) 



-2n 



1 

'20 



ca cp (2c — 1) 



■ + 



2/i2 2/i 4/i2 



+CU I cqe 



V 

2^ 



z/(z/ + 8cg) ?7 



32/i 4^2 



+ e" 



2V2cV-3r/ z/2 -8(2cp + c2g2) _4c2a2 



4^2 



+ 



32/i 



c^g^ 2 - u / 1 2 / 2 N 2 

^ + ^(cga + -.^ + (c^-c),^ 



(4cV + 3cV-i^')^^^- (1-33) 



Theorem 1.2. ^ As n ^ 00 



c{q{s) + n(s))e-'^(^) - ^(l - 6"^^) 



Ms) 



n 3 + 



(1.34) 



uniformly for hounded s . 

1.1 Statement of our results 

If we set 

/oo 
qn{x) dx and h{t) 

then for the Gaussian Orthogonal Ensemble, 
Theorem 1.3. For n even, 



Pn{x) dx 



1 2 m o Kt) 



2 '''ait) "V2a(t)^'' 
1 , „2 Kt) \ ^^gj^ ^2a{t)h{t) 



sinh ^2a{t)h{t) 



2 ''ait) 

and for the Gaussian Simplectic Ensemble, 
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Theorem 1.4. For n odd, 

F^4t/V2Y = F„,2(t)^ (l + cosh ^2a{t)b{t)^ 

Theorem (11 .41) leads to the following Edgeworth Expansion of F„ 4 
Theorem 1.5. Then as n — j- 00 

FnA^? = Ms) ■ {cosh\^) + |[no(l + cosh(/x)) - gsinh(/i)]n-3 



u , , , r, o , , , cosh(u) — 1 , , 

sinh(/i) + c^q^ cosh(/i) + Ec,2 + (1.35) 



.2x/2/x '"^^ '"^^ 10 

•\/2(?7 — A/2c^gno) sinh(/i) 
uniformly for bounded s . 



In §2 we complete the derivation of fll.24p as outlined in [22]. In §3 Theorems 11.31 
and 11.41 are derived. In §4 we justify Theorem 11.51 followed in §5 by a brief discussion 
on the rate of convergence of these distributions as a function of the fine-tuning 
constant c for large n. 

2 Derivation of 4 

We treat here the case n odd. (Most of the algebra used here can be found in [TO] . 
[H], [12], [S], and [27].) We note here that 



FnAit/V2y = det(J - K„,4), (2.1) 

with 

Kn,2 + e Sip Kn,2 + eip®ij J ^ ^ ' 

We set J = (t, 00), and x represents the multiplication by the function Xj{^)- This 
notation allows us to think of the operator with kernel Kn^A as an operator acting on 
M instead of acting on {t 00). 

We will like to also note that this section is identical to the derivation made in [22] . 
we only provide here details intentionally left out by the authors of [22] (we suppose to 
shorten the length of the paper, even if they provided the necessary steps to complete 
the derivation.) 

Using the following commutators, 

[Kn,2, D]=(p^1p + 1p^ip, [e, Kn,2\ = ^ Elp - Elp ^ Eip (2.3) 
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{ip and (f appear as a consequence of the Christophel Darboux formula applied to 
Kn,2i) we have 

Kn,2 + i'®£'^ = D eKn,2 + D eip ^ Eip = D{eKn,2 + si' ® sv) = D{Kn,2£ - ® 

Kn,2 D -ij®^ = D Kn,2 + ^ ® = DKn,2 + De^ ® = D {Kn,2 + e<^ ® tp) 

and 

£ Kn,2 + Sip ® e^p = Kn^2 E - E If ®) Elp 

as De = I. Our kernel is now 

' V Kn,2E-Eip®E1p Kn,2 + E^®i) I ^ ' 



(2.5) 



XD 0\ / {Kn,2E - Eip®Ei))x {Kn,2 + Eif ® i))x 
X ) \ {Kn,2 E - E If (g) Elp) X {Kn,2 + Eif> ® ^p) X 

Since A'n,4 is of the form AB, we can use the fact that det(/ — AB) = det(/ — BA) 
and deduce that the Fredholm determinant of Kn,i is unchanged if instead we take 
2Kn 4 to be 



{Kn,2E-E^®Ei))x {Kn,2 + E^®^)x\ ( XD 
{Kn,2 E - E if (g) Elp) X {Kn,2 + Eip ® 1p) X J \ X 



(2.6) 



{Kn,2E - E^®E%p)xD {Kn,2 + E^ ® ^)X \ /2 7) 

{Kn^2E-E^®E^)xD {Kn,2 + Eif ® ij) X I ' 



Thus 

det(/-K„,4) = det 



l{Kn,2E-Eip®Etp)xD -l{Kr,^2 + E^ (S) ij)x 
-\{Kn^2E - Eif® Ei))xD I - \{Kn,2 + ELp®i))x 

Performing row and column operations on the matrixo does not change the Fredholm 
determinant. We first subtract row 1 from row 2, next we add column 2 to column 1 
to have the following matrix 

I - \{Kn,2E-Eip®Etp)xD- l{Kn,2 + E(p ® tp)x -U^<n,2 + E(p®1p)x 

/ 

(2.9) 

We therefore have, 

det(/-K„,4) = det(^I -^{K\,2E-E^®Eij)xD-^{Kn,2 + E^®ij)x^ (2.10) 

= det(^I - Kn,2X-\iK^,2E-E^®E^)xD + ^{Kr,,2-E^®ij)X^ (2.11) 

= det f / - Kn,2X - ^ {Kn,2 E-E(f®E^)xD + ^Kn,2X + ® i^X - £f ® ^^X 



^This does not change the determinant, for more details see [22] 
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= det ^/ - Kn,2X - \^n,2 {sxD - x) - ^e<^ ® ip{exD - x) - ® 
= det - Kn^2X - \^n,2 {sxD -sDx)- ® ^(^X^ - £ D x) - e^p ipx 

= det ^/ - i^„,2X - \^n,2 e[x D] - ]^eip ® ipe[x D] - Sip ® i^x 

We used the fact that e D = I to write x = ^ -D x and the fact that e is antisymmetric 
to have 

Sif ® eil) X D = Sip ® ipe^ xD = —e(p ^ipexD. 

Next we factor out / — Kn,2 and note that (/ — Kn,2)~^ = I + Rn,2, where Rn,2 was 
defined as the resolvent of Kn,2, and (J — Kn,2)~^£^ = Qn,e- We are interested on the 
determinant of the following operator 

{I - Kn,2x){l -l{Kn,2 + Rn,2Kn,2)e[x D] - ^Q^,, ® e[x D]-Qn,e®^x)- (2.12) 

The determinant of the first factor, det(J — ii'„^2X) = ^n,2 is a familiar object that 
was studied in our work in GUE„ see [2] or [3] . Attention will be given to the second 
factor of fl2.12p . We will represent this factor in the form (J — X]j=i ® f^j) and use 
the well known formula det (J — Yl^=i ® l^j) = det {5ij — (a^, ^ to expand 

the Fredholm determinant. First we need to find a representation of e[x-D] as a finite 
rank operator. To this end we introduce in this section the following notation, 

ek{x) = e{x-ak), Rk{x) = Rn,2{x,ak), 6k{x) = 6{x-ak), ai = t, and a2 = oo. 

With the new notation J = [t, oo) = (ai, 02), and the commutator 

[X D] = -61 ^61 + 62® 62, 

gives 

e[x D] = -El (g) 61 + £2 ® <^2- 
With this representation, we have 

{Kn,2 + Rn,2 Kn,2)e [X D] = Yl + ^",2 Kn,2)ek ® 4 

k=l,2 

and 

k=l,2 k=l,2 

In the last equation we used the formula (a ® /9) ■ (7 ® 5) = {f3,'j)a ® S. If we 
substitute this in the second factor of (12.121) . it becomes 



I-7^Y1 (-l)''(^n:2 + Rn,2Kn,2)£k ® ^fc - - ^ {-if {^J, €k)Qn,e ® Sk - Qn,e ® X^- 

(2.13) 



k=l,2 k=l,2 
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We have 



1 1 

£2{x) = Eooix) = -sgn(x - oo) = --, and R2 = Rn,2{x, 00) = 0. 



If we substitute these value in fl2.13p . it then becomes 

[{Kn,2+Rn,2 Kn,2)et+ {ip , et)Qn,s] ^St + -^ 



I-Qn,e®X^-\ [{Kn,2+Rn,2 Kn,2)et+ {tp , et)Qn,s] ^St + ][{Kn,2+Rn,2 i^n,,2)l+ l)Qn,e] ®5o 



(2.14) 

This operator is of the desired form 

3 

/ - ^ ttfc (g) /3k and det(/ - ^ ® = det{5ij - (a^, /^i))^^^-^! 3 

fc=l,2,3 j=l 

with 

(^1 Qn,£i 

aa = ^ [{Kn,2 + Rn,2 Kn,2) + (^A, ^i) QnJ , 
as = -i[(i^n,2 + Rn,2Kn,2)l + (V^, l)Qn,.] 

and 

/3i = x^, = (5t and /^s = Joo- (2.15) 
As pointed out in [21], the contribution of fi^ is zero 

("1 , /^s) = ("2 , /^s) = ("3 , Pz) = 0, 

thus the determinant reduces to the contribution of ai, 02, /3i and /32- The corre- 
sponding inner product are; 

(ai, /3i) = (ai, /^a) = + c<^ = qn,e as c<^ = for n odd (2.16) 

where we have 

1 [°° 

c^p = eip{oo) = - (/?(x) dx = as 99 is odd (2-17) 

2 j-00 

we also note that 

1 f°° 

c^ = eip{oo) = - / ip{x)dx, (2-18) 

2 J-CX) 

and for n odd 

c, = (vr (n - i))i/42-3/4-(n-i)/2 ((^-l)0^^^ ^ (2.19) 



and 



(tta, /3i) = -^(K,4 - ^4 + a^Vn,e), («2, /32) = -^(^n,4 + a4g„,£). (2.20) 
with 04 = (■?/; , Et). The determinant of fl2.13p is therefore 

(1 - VnAt)) (1 + ^^n,4(t)) + ^gn,e(t) K,4(t). (2.21) 

We will make use of large n expansion of the functions appearing fl2.2ip (see [5]) 
combine with the large n expansion of F„^2 to give an edgeworth type expansion of 
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3 Finite n Expression of i and 4 

In [5] we derived the following results; if we recall that 

(^i't) = / Qn{x)dx and b{t) = / Pn{x) 
Jt Jt 



dx 



then for the GSE„, 

in,eit) = 1 - ^[1 + cosh ^2a{t)h{t)] (3.1^ 

and 



<ln,eit) = -J^sinh v/2a(t)6(t) (3.2) 



for the epsilon independent functions, and the epsilon dependent functions are 

(3.3) 
and 



VnAt) = -c^^[l + cosh ^2a{t)b{t)] - J sinh ^2a{t)b{t), 



TlnAt) = -c^,y ^ sinh ^/2a{t)bit) + cosh ^/2a{t)b{t) - 1. (3.4) 
If we recall that 

-Fn,2(t) = det(/ - i^„,2) = exp ^-2 ^ (x - t)g„(x)p„(x) dx^ , (3.5) 

then have 
Theorem 3.1. For n odd, 



F„,4(t/V2)' = i"n,2(t)^ (l + COsh 



Or what is the same as 



FnA{t/V2) = cosh ySMSexp (x - t)g„(x) p„(x) dx^ (3.6) 

The surprising fact here is that drops out and the formula is as simple as pos- 
sible. We couldn't wish for a better result. 

In a similar way we can also use the following representations 

Fn,l{tf = F„,2(t) |(1 - Vn,e){l - ^7^„,l) - ^{qn,e - C^)^n,l| • 

together with the following representations (see [5] for a derivation of these results) , 
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:it) = 1 - i[l + cosh x/2a(t)6(t)] + c^J ^ sinh ^2a{t)b{t) (3.7) 



and 



qnAt) = -^/^si"^^ ^/2a{t)b{t) + c^cosh ^2a{t)h{t). (3.8) 



We also have 



Vn,i{t) = c^^[cosh ^2a{t)b{t) - 1] - y ^ sinh ^2a{t)b{t) 



(3.9) 



and 



7^„,l(^) = 1 - 2c<^ sinh v/2a(t)6(t) - cosh ^/2^^(t)b(t) , 



(3.10) 



to have for the GOE„, 
Theorem 3.2. For n even, 



Fl,(t) = Fn,2{t) { \ - 4^ - 2 J sinh ^2a{t)b{t) 



1 , ^2 Ki) 



2 '^a(t) 



^ cosh v/2a(t)6(t)| (3.11) 



Unhke the GSE„ case, the finite n GOE„ does not simphfy as well, but it is still 
a very applicable formula. 

Theorems 13.111 and 13.61 are the main results in this section, they provide a general 
n formula for the probability distribution function of the largest eigenvalue. In the 
next section, we will find a large n expansion of these two functions. 

4 Edgeworth Expansion for GSE^ 

The previous section tells us that these three probabilities distribution functions 
Fn,i, -F„2 and F„4 are all functions of g„ and thus to find a large n expan- 
sion of these functions, we need large n expansion for these functions g„ and 
Upon substituting these expansions into the sought after functions (13.111) and (13. 6p . 
we recover an Edgeworth type expansion of these probability distribution functions. 
This work was carried out in [21 [3] for the GUE„, and in |1] for the GOE„. In this 
paper we follow the same technique to find the corresponding expansion for GSE„. 
This presentation is different from the previous results since it gives a closed form for 
Fn^i and Fn^A that can be used to compute these expansions. 
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We quickly see that after substitution of the large n expansion of g„ in a{t) and 
the large n expansion of p„ in the result is 



i^„,i(r(.)) = F„,2(r(s)) ■ |e-^(^) + (^^(e'^W - 1) + cg(.)e-'^Wj n't 

(e-'^(^) (z/(s)' ((-1 + e^(^)) (-3 + 8c + e^^^^) + 2/i(s) (l - 4c + 4cV(s))) 



16/i(s) 

8cz/(s) (-1 + e''^') + + 2c/i(s))) / rfx + 



8/i(s) (^20c (-1 + e'^^'^)) q[x\v[x\ dx - j qi[x\ dxj + /i(s) ((-3 + 20c2) 

X) poo poo poo poo 

p[x]u[x]dx—3 / (ia;— 2 / i;[a;]pi[x] rfx— 2 / p2[x]dx—3 / qi[x]dx+ 



r'OO /poo 

c^ I ( / (ix)^ — 20 I 2 / dx — 3 I q[x]v[x] dx + I qi[x] dx ] ] + 



2 



M[x]g2H'^a^+ / c/x + J q[x]u2[x]dx— j p[x]vi[x]dx] ]]]n ^ 

+0{-) 
n 



^ ' (4.1) 
Using the Edgeworth Expansion of F„ 2 given in Theorer dl.il leads to Theorem ll.2[ 
For the GSE we have the following, 

Fr^Ai^ f = FnAr{s)) " {cosh^l^) " ^cg(.) sinh/.(.) 



1 



f 00 

4c^ cosh[/i(s)]/i(s)(z/(s) — / 



16/i(s) 

-z/(s)^ + 4/i(s) ( (3 - 20c^) / p[x]M[x](ix + 3 / 



00 poo poo / poo 



2 / (ix + 2 / p2[x]dx + 3 / qi[x]dx + 20c (2 / g[x]M[a;] dx— 

Js J s J s \ J s 

poo poo \ / poo poo 

3 / c/x + / gi[x](ix) — 21 / ■u[x]g2[a;] c^a; + / gi[x]Mi[x] 



00 poo 



q[x]u2[x] dx — J p[x]vi[x] dx j j j sinh.[fi{s)] j n 3 + 0(n ^) 

= F„,2(r(s)) ■ |cosh2(^) - ^cq{s) smh fi{s) n-l +g{s)n--^ +0{n-^)^ 
And if we substitute the Theorer dl.il we recover Theorer dl.SI 

12 



(4.2) 



5 Conclusion 



The major finding about this studies of the probabihty distribution function of the 
largest eigenvalues of the classical Gaussian Random Matrix Theory Ensemble, is 
their dependence on just two quantities, g„ and p„. Indeed, comparing the represen- 
tation for these three ensembles, we see that a complete statistical study would only 
require the knowledge of of these two functions. The simplicity of representation fl3.5p 
and fl3.6p allows us to have a simple form for the probability density function of the 
largest eigenvalues, and the somewhat not so simple representation (13.111) will do the 
similar work for the GOE„. 

The second finding is the similarity between the unitary and the Symplectic ensem- 
bles, the probability distribution functions are very simple to represent, the fine- 
tuning constant c does the work, as in these two cases, we can eliminate the first 
correction term to the Tracy- Widom limit by setting c to zero. But there are some 
numerical evidences showing that there is a way to fine-tune the result for the GOE„ 
to speed the convergence rate to the Tracy- Widom distribution. This gives us hope 
to find the right way to capture this phenomenon for the Orthogonal Ensemble. 
It would also be convenient to find a simplify version of the functions f{s) appearing 
in (14.11) and the function g{s) appearing in (14.21) . 
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